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00 ■ Abstract 

We construct the Gibbs state for z/- dimensional quantum crystal with site dis- 
placements from ]R , d > 1, and with a one-site non-polynomial double- well 
potential, which has harmonic asymptotic growth at infinity. We prove the 
uniqueness of the corresponding Euclidean Gibbs measure (EGM) in the light- 
mass regime for the crystal particles. The corresponding state is constructed 
via a cluster expansion technique for an arbitrary temperature T > 0. We 
show that for all T > the Gibbs state (correlation functions) is analytic 
with respect to external field conjugated to displacements provided that the 
mass of particles m is less than a certain value m* > 0. The high temperature 
regime is also discussed. 
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1 Introduction 

It is generally excepted that for investigation of different physical phenom- 
ena in quantum crystals one can consider an infinite system of interacting 
anharmonic oscillators, which are situated in the sites of //-dimensional lat- 
tice Zi v . The heuristic Hamiltonian of such a system (in the case of 2-body 
interaction) has the following form: 

H = -4: E A ^' + E w fa) + \ E Qb*> %) ~ E H> 

jez v je% v i,je% v ,i^j 3^% v (1-1) 

where m is the mass of particles, the operator Aj corresponds the kinetic 
energy of the system and, in fact, is d- dimensional Laplace operator in the 
one-particle Hilbert space L 2 (]R ,dq), where dq is the Lebesgue measure on 
M , qj G M for j G 5Z U is displacement of a particle from its position in 
the site j G 7Z V . For general case d < v > 1. But the most of results, 
which where obtained earlier are for d — 1. The particles are confined near 
their sites by potential W(qj). A harmonic one-site potential Wharm{Qj) = 
\aq 2 j, a > together with a harmonic two-particle interaction Q(-, •) define a 
well-known harmonic crystal model (1.1). To produce a model describing a 
(ferroelectric) structural phase transition one usually takes for W(-) a double- 
well anharmonic potential, keeping Q(-,-) harmonic (see e.g. [1,19,24]). For 
example, 

W(q) = 1 -b(q 2 ) 2 -±aq 2 , a,b G R\, 

or semibounded from below polynomials of higher degree with more than 
two equal minima. Here q 2 = q ■ q is the scalar square of displacement vector 
q G M . The proof of existence of phase transition in such kind of systems 
for harmonic interaction Q{qi,qj) were obtained in [12,15,26,34,40,47,53] for 
the case 



Q(Qi, Qj) = QijiVi - Qj) 2 , (1-2) 

where Q = (Qij)ije^ v i s a matrix of non-negative harmonic-force constants. 
The term h ■ qj = X]»=i h a qf break the symmetry in the direction of the 
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external field h. The non-uniqueness of Gibbs states for h = is proved 
when the mass of the particles m is sufficiently large and the temperature T 
is sufficiently low, or (3 = (kT)" 1 is sufficiently large. 

From the other side, another interesting phenomenon, the suppression of 
the long-range order by strong quantum fluctuations in such systems was 
experimentally observed (see, e.g. [55]) and was discussing long time ago 
from the physical point of view, see [51], or the books [1,19]. A rigorous 
study of this phenomenon was given in the series of papers. For the d = 1 
displacement and all v > 1 the suppression of the order parameter for all 
temperatures including T = (as soon as the particle mass is smaller than 
some threshold: m < m*) was proved in [58]. Later it was shown in [3] 
for d = 1 and in [37-39] for d G IV, that not only the long-range order but 
also any critical behaviour of displacements of particles from the equilibrium 
positions are suppressed at all temperatures T > if the model is "strongly 
quantum" , i.e., if the mass m is small enough. But a more important problem 
is to prove the uniqueness of the Gibbs state for such systems for sufficiently 
small mass of particles. It would imply the absence of all critical anomalies 
in the regime of "strong quant umness" . The first step in this direction was 
made in [44] for d — 1, where for small values of m the convergence of 
corresponding cluster expansions was proved for all values of temperature 
T > 0. See also [8], where the cluster expansion is constructed by another 
and simpler technique but for small m, which depends on the positive T > 0. 
The result of [44] has been extended to the multi-dimensional spin in [27]. 
This in particular implies that the order parameter is always zero if m < m* 
for any values of T > [43]. It gives, also, uniqueness of the EGM (see 
[2,31,33] and [7] for detail description of the Euclidean approach) at least for 
so-called "compact" boundary conditions for displacements. In particular, 
this is true if displacements are reduced to a compact (so-called compact 
spins). This resumes the results of [44] and [43] about the suppression of 
the displacement phase transition by the large quantum fluctuations in small 
masse domain m < m*. For models with the compact spin the uniqueness 
was proved in [9]. Then in [4] the uniqueness of EGM was proved for d = 1 
and for every fixed inverse temperature (3 < oo if the mass is less than some 
threshold: m < m*(f3). Finally, the dependence of the mass threshold on the 
temperature was eliminated in recent papers [5], [6] for the set of tempered 
EGM. 



Remark 1.1 These results were obtained for the case of polynomial anhar- 
monic one-site potential of the form: 

s 

W(q)=v(q 2 ) = Y,a2 P (q 2 ) P , a 2s > 0, s > 1, s6ff. (1.3) 

P =i 

with additional assumption on W(q) that v(-) is convex on M\. For the proof 
of these results the fact of smallness of the parameter ?77,( s-1 )/( s+1 ) for small 
m, i.e. s > 1, is crucial. Recall that in this case, after an appropriate unitary 
transformation of the Hamiltonian, the parameter m^^ 1 ^* '( s+1 *> stands in front 
of the harmonic interaction term Q. So the light-mass regime corresponds 
to the one-site energy domination and the weak interaction of anharmonic 
oscillators defined by W , that implies the uniqueness of the EGM. 

On the other hand, let us define s > by 

s = inf{a : limsup v(q 2 )/(q 2 ) a < oo}. (1.4) 

If s < 1, then in the light-mass domain (" strong quantumness") the param- 
eter 777,( s_1 )/( s+1 ) is not small, i.e. the approach developed for s > 1 does not 
work. So, the rate of the growth of anharmonic potential is an important 
parameter distinguishing different regimes, see discussion in [43] and [58]. 

Remark 1.2 If s < 1, then in the light-mass limit the harmonic interaction 
term dominates the one-site anharmonic potential, i.e. the system has to 
behave close to a Debye harmonic crystal (not proved). In the marginal case 
s = 1 the both harmonic interaction and harmonic part of W dominate the 
anharmonicity. From the physical point of view the large-distant behaviour of 
the one-site anharmonic potential is not well-justified. So, a priori all above 
scenario are possible. 

For example, the paper [54] proposed a model with the one-site non-polynomial 
anharmonic potential of the form: 

W(q) = -aq 2 + be' 5q2 / 2 , a,b,8 e R\. (1.5) 

corresponding to s = 1. It is easy to see that for a < bS the minimum of 
the potential (1.5) is not unique and the system may undergo a displacement 



phase transition. For the case of so-called self- consistent phonon approxima- 
tion [1], or for the displacement quantum spherical model it was proved in 
[48] for (1.3) and in [54,58] for (1.5). In general case this follows from [15], 
see also [47]. On the other hand, in [48] and [54,58] it is shown that there is 
a non-zero threshold mass m* such that the phase transition in the displace- 
ment quantum spherical model is suppressed by the "strong quantumness" 
as soon asm< m*. 

The aim of the present paper is study the case of the quadratic growth, 
s = 1, of the double-well anharmonic one-site potential W (1.5). We prove 
that (similar to the case s > 1) there is a non-zero mass m* such that 
in the light-mass domain: m < m* the critical behaviour of the quantum 
anharmonic system (1.1) for h = is suppressed for arbitrary dimension 
of particles displacements d > 1 and v > 3. More precisely, we prove the 
uniqueness of the corresponding Gibbs state. Our result is based on the 
proof of convergence of the cluster expansions for the quantum Euclidean 
Gibbs Measure in the thermodynamic limit for arbitrary temperature T > 
and sufficiently small mass m < in* of particles (uniformly in temperature). 
Notice that since for s = 1 the interaction in the light-mass limit m < m* 
is not negligible (Remark 1.2), our strategy is different to that for s > 1, cf 
[44] . In that sense our result supplements the recent proof of the Gibbs state 
uniqueness for s > 1, see [5,6]. 

Now there are few other remarks in order to outline our strategy of the 
proof the uniqueness for the case s = 1 . 

The anharmonic perturbation in (1.5) can be re-written in the following 
form 

V(q) = be- Sq2 / 2 = b I dfi(a)e tV ~ Sa - q , (1.6) 



where d/j,(a) = (2ir)~ d / 2 e~ a l 2 dct is the Gaussian measure on M, . So anhar- 
monic term (1.6) is a gentle perturbation of the quantum harmonic system 
in the sense of Albeverio-H0egh-Krohn [2]. Unfortunately, our situation is 
similar to [2] only at the first glance. In fact, there are some essential differ- 
ences. In [2] the bounded measure d\i has to have a compact support. Under 
this condition the authors managed to apply the Kirkwood-Salsburg Equation 
(KSE) technique to prove the existence of the corresponding correlation func- 
tions in thermodynamic limit as well as their analyticity in coupling constant 
A (in our case A = b) in some domain |A| < Aq for all temperatures. 



We can not apply this technique, since of the Gaussian measure in (1.6) 
has infinite support in ]R d . By consequence it is impossible to find a Banach 
space E% for correlation functions and the corresponding KSE. As it was 
indicated in [2], Remark 2, for extension of the results to the measure with 
unbounded support one needs a modification of the Banach space. But even 
that is not essential, since the main difficulty for s — 1 is connected with the 
parameter of analyticity. Instead of A we have as the small parameter the 
particle mass to, which is involved in the Hamiltonian via the kinetic-energy 
term, see (1.1). On the other hand it is impossible to expect the analyticity 
in domain \m\ < to*, because we must always consider Rem > to define 
the system (1.1). 

We overcome these difficulties by a cluster expansions technique. The 
choice of the cluster expansions is defined by the dependence of the param- 
eters in the Hamiltonian on the mass to after an appropriate rescaling, see 
(2.6). In the case of anharmonic potential with harmonic asymptotics s = 1 
the dependence on the mass to remains only in the rescaled temperature and 
in the rescaled anharmonic potential (1.5). Therefore, the expansions should 
be arranged with respect to (w.r.t.) the anharmonic part of the one-site 
potential. (Recall that for s > 1 the corresponding cluster expansions are 
constructed w.r.t. two-body harmonic interaction.) This gives a possibility 
to construct the limiting Gibbs state and to prove its analyticity in some 
small circle \m — to*/2| < m*/2 of the complex values of to around the point 
to*/2>0. 

To prove the uniqueness of the corresponding EGM we profit the Theorem 
3.1 of [6] and we use the cluster property of the EGM moments, which is a 
consequence of the cluster expansion convergence established in the present 
paper. Similar to [43] we prove the analyticity of the state w.r.t. the external 
field h. As a consequence this immediately implies that the displacement 
order parameter (in the region of convergence of the cluster expansions) is 
equal to zero for all temperatures T > as soon asm< to*. 

A brief contents of this paper is the following. In Section 2 we describe our 
system (1.5) in more details and we formulate the main results. In Section 3 
we construct the cluster expansions w.r.t. the anharmonic part of the one-site 
potential, and we give the proof of main theorems. In Section 4 we provide 
all necessary estimates. Section 5 is devoted to the limiting Gibbs state 
analyticity w.r.t. to the external field, and in Section 6 we briefly discuss 
the proof of the convergence of cluster expansions and the uniqueness in the 

5 



case of high temperatures. 



2 Description of the system and main results 

We consider a u(> 3)-dimensional quantum crystal model with d(> 1)- 
component site displacements (or d- dimensional continuous infinite-spin vari- 
ables) on the cubic lattice 7Z, V . With each site j G /Z u we associate a one- 
particle physical Hilbert space L 2 (]R d , dq), where dq is the Lebesgue measure 
on R d . Then 

H A = 6^L 2 (R d ,dqj) ~ L 2 (M d \ A \dq A ), 



jeA 

dqA = J\dqj, 

jeA 

is the Hilbert space related to some given bounded set A C Z v , with number 
of sites |A| < oo. 

For a finite set A C Z v we consider in 7i A the Hamiltonian: 

H A = < + V A . (2.1) 

Here the harmonic part of the Hamiltonian has the form: 



jeA jeA (i,j) C A (2.2) 



where m is the mass of particles (we put h = 1), Aj is the Laplace operator in 
L 2 (IR d , dqj), the sum is extended over all pairs (i,j) C A for which \\i— j\\ = 1 
and J > 0. We consider a nearest-neighbors interaction only for simplicity. 
The constant C\ we will be chosen later. 

We put a one-site potential W(q) in (1.1) to be equal to (1.5) for q 2 = q-q, 
i.e., for q £ 1R . So the anharmonic part of the one-site potential in (2.1) has 
the form 

V(q) = b e~ Sq2/2 , b>0,5>0, qeR d . (2.3) 

To prove the existence of the Gibbs states in the light-mass regime for 
small temperatures, including T = 0, we need a supplementary argument 
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based on the cluster expansions. Below we do this first for the periodic 
boundary conditions (p.b.c), then we complete our construction by the case 
of the empty boundary conditions The last we use later to pass to tempered 
non-empty boundary conditions Notice that the cluster expansions for p.b.c. 
prove also the clustering property of the corresponding EGM for m < m* and 
for any T > 0. Therefore, simultaneously with existence we get uniqueness 
of translation invariant EGM for all temperature, including T = 0, as soon 
as the particle mass is light enough: m < m*. So, first we assume in (2.1), 
(2.2) the p.b.c: 

A = {j'Gr| -NJ2 < jM < 7V M /2, /i = 1, ..., v- NJ2 e IN}. 

Thus, |A| =N 1 ---N V and 

^),,i(")+^,.,iM) = q t^ a = *' -' d ' ^ = 1 ' -' v - ( 2 - 4 ) 

Then, to complete the proof of the uniqueness of the EGM, in the next 
section we consider also a general boundary conditions generated by some 
fixed configuration oo from the set frL, which includes so-called tempered 
configurations, see [6,13,14] and [46]. 

Let £(7^a) be algebra of bounded operators on Ha- We consider the 
temperature quantum Gibbs state on C{7i\) defined by 

Pa[A) ~ z p (A) ' p ~ k B r (2 - 5) 

where Z P (A) = Tre"^* and A e C(H A ). 

Note, that the algebras C{Tl\) (for different A) are isometrically embed- 
ded into C{7i\>) for A C A'. Using this fact we can define the inductive 
limit 

2l = lim C(Ha), 

which a normed *- algebra of local observables. The closure of this algebra 
in the operator norm forms the algebra of bounded quasi-local observables 

21 = %. 

The main result of the present paper is the following theorem. 
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Theorem 2.1 For the system of quantum particles with Hamiltonian (2.1)- 
(2.3) there is some m* > ; such that for any < m < m* the limiting Gibbs 
state 

lim pl(A) = pP(A), Ae%, 

- exists for all T > 0; 

- it is a unique translation invariant quantum Gibbs state (corresponding to 
EGM) on the algebra 2lo, which can be by continuity extended to the algebra 
21. 

Notice that by virtue of results of the paper [32] (see also [7,33,44]) it is suffi- 
cient to prove this theorem for some sub-algebra 2T C 21 of local operators, 
which we describe below. The main technical tools are the Feynman-Kac 
formula and the representation of the Gibbs states (2.5) on this algebra by 
functional integrals with respect to some measure //(•) which we are going to 
construct. 

Since we are going to consider domain of small masses, the following 
standard light-mass rescaling of variables (appeared for the first time in this 
context in [58]) is relevant: 

qj = aXj, a = m~4. (2.6) 

This change of variables induces the unitary map : 

U : Ha — > Ha 
(Uf){x A ) = a lAl/2 f((ax) A ), (ax) A = {axj, j E A}. 

One can check that 

UHaU- 1 = m-^H A , 

where 

H A = H£ + V A , (2.7) 

Ho = -\Y.^ + \ a H x ) + J E (^-XjY + CL 

jeA jeA (i,j)cA (2.8) 

Ci = Caw 1 ' 2 , (2.9) 



V A = b m J2e~ Smx2j/2 , b m = bm l '\ 5 m = 5m' 1 / 2 . 

jeA (2.10) 

We also define a rescaled inverse temperature 

P = p m -^ (2.11) 

which yields 

U{3H K U- 1 = J3H A . 

Then the Gibbs state (2.5) takes the form: 

A{A) = p{(A), A = uau-\ 
where 

p{(.) = Z7\A)Tr{.eM-PH A )} . (2.12) 

To represent the state (2.12) in the form of the functional integral we 
rewrite the operator H A as follows 

^o A = E^o a) ' A , (2-13) 

a=l 



Hi a)A = ~\^ + \(^\ Bixf) - l -TrB A , (2.14) 



where 



A a =EtW B 2 A = -JA A + at A , (2.15) 

jeA ° x j 



x a = ( x j )jeA £ JR , and B\ is a periodic |A| x |A| matrix, defined by the 
elements: 



2v, for j = k, 

A A ).;/,= <; -1, for |j - fc| = 1, (2.16) 

0, otherwise , 



9 



and 

Wjk — Sjk = 8jWkW ' ' ' ^jMfeM- 

Now we fix the choice of the constant C A in (2. 2), (2. 9) by 



C A = ^TrfiA. (2.17) 

2\/m 



By virtue of p.b.c. the lattice operator B A is cyclic (see (2.16) and [30, 
Sec. 16.4]), so we can calculate its eigenvalues explicitly: 



\ k = y/I(k), £ (k) = a + 4J^sin 2 (A;^/2), fe e A*, 

m=i (2.18) 

where the set A* is dual to A, i.e. 

2ir 
A* = {keM"\k = {k^}; =1 , k^ = —n^\n^ = 0, ±1, ..., ±N^ - 1/2, N^/2}. 



'ii 



(2.19) 



Now the Feynman-Kac formula (see e.g. [50], [52]) gives 
Tre -/3H A = f dx A (e~^ A 

JjRdlAI V 






^a / <;JN^ J »^ A(T) » (2.20) 

7iR d l A l 

with 






£>(x A ) = \{x K , B 2 A x A ) - ^TrB A + V(x A ), (2.21) 



where we put V(x A ) = V A , see (2.10), and 

wi, XA (dcv A )^]JW^ Xj (d 

jeA 



Uj), 
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here W£.. x .(duj) is the conditional Wiener measures on the space of contin- 
uous periodic trajectories 

n aA = K(-)K : 5 a - iR d|A| }, 5, := [0,/3], 



"13, A — l^A\ j\^a ■ ^p ' -"«- j, ^ 

2g, iZ^ := o^ 



w A (0 = K-(-), J e A\uj e fy, fy := C(S* - iR d )}. 



Now, following [2] the unperturbed measure in the right-hand side of (2.20) 
(i.e. V A = 0) can be rewritten with the help of a Gaussian integral over the 
Hilbert space: 

% := K e (% A , S M ) I |M| 2 < oo}. (2.22) 

Here Ea a is the standard a-algebra of fig A -subsets generated by Borel cylin- 
der subsets and norm is defined by 



,,,, / " /3 |-,gL>a(t) du\(r) 
\ u a\\ = I clt 







(-^,-^) + MT),a; A (T)) 



where (•, •) is the scalar product in ^(M ' '). Then for any real continuous 
function F defined on Ti^ A with this unperturbed measure one gets: 

dXK I Wl. x (doj A )e-yo^ A (r),B^ A (r)) + §TrB AFM = 



. . , x A ;x A 

iR dA Jn- 



= C [ d//£(0 A I P)F{4>a) = Co [ rf/i^(0 A )F(0 A ), 

Jh 0,a J K0,a (2.23) 

where Co is a normalization constant and dfiQ (0 A | p) is the Gaussian measure 
on Tig A which corresponds to p.b.c.. We skip the index p in the following. 
By (2.20) and (2.23) it is clear that 

C = Z?(A) = Tre"^o A . ( 2 .24) 

Hence the right-hand side of (2.23) is the expectation: 

(F)^ = I rf/4(0 A )F(0 A ) (2.25) 



W /3,A 
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with respect to the Gaussian measure on the Hilbert space of continuous 
periodic functions generated by the corresponding Gaussian process with the 



zero mean and the covariance Gq. 



G^ir -r')= d^{Mf\T)<t>t'\r') = S«a'Gk dk (T - r'). 



W /3,A 



Here a, a' E {1,2,..., d}, and G^.^t — t') is the Green function of the oper- 
ator (-£2) 1a + B\ on L 2 (A <g> S^ R d ) with p.b.c. More explicitly [50]: 

r A , j^ _ 1 v-^ A3~k)-i l V" cos(27r(r - r')n/(5) 

l A '^ /^ (4vr%V^ 2 )+^) (2.26) 

or taking the sum over n, we get (see e.g. [50]) 



^'-o-iliE^ 



e (/3-|r-r'|)V^(0 + e k-rV £ « 



l A l^ 2 v ^(0(e^v / ^)-l) ' (2.27) 

It is important to note that Gq depends on the mass m only via rescaled 
temperature (3 = ftm^ 1 ' 2 . 

Remark 2.1 From (2.26)- (2.27) one gets that for A / Zl> v the covari- 
ance Gq. - fc (r — t') converges pointwise to a translation invariant function 
Gq.\j-u{t — t') which is also a bounded continuous function and in fact it 
is the kernel of a positive self-adjoint operator in L 2 (5Z U <8> Sz;]R ). From 
this convergence it follows that the corresponding process (measure) converges 
weakly to a limit process, which is a homogeneous Gaussian process on *Z V 
and periodic on Ss with values in lR d . 

By virtue of (2.20) and (2.23) it is clear that 

Z A = Tr A e"^ A = Z° A [ d^{<P)e~ ^ dr ^^\ (2.28) 
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Now we can follow the line of reasoning of [2,7,32]. For every bounded func- 
tion A(xa) on lR d > A > we consider a bounded operator A on Ha defined as 
multiplication on a bounded function: 

(A f)(x A ) = A(x A )f(x A ), (2.29) 

and for any t > we consider the operator 

A t = e- tAA A e ttiA . (2.30) 

Then for every set of bounded functions A^ ' , . . . , A^ 1 ' and for increasing 
sequence of the "time" moments = to < t\ < . . . < t n < t n+ \ = (3, we 
define the operator 

n 
A _ TT^W — ^4(°) e -ti^A^4( 1 ) e -(*2-ti)^A e -(t n -t n - 1 )H Kj ^n) e t n H K 

1=0 ' (2.31) 

Then by (2.12) one gets: 

pP(A) = 1 Tr(,l (0) c~ flJ?A /l (1) c~ (f2 ~ fl) ^ A . ..e-fc-Vi^A^Wg-lHt.)^^ 
Z ^ ' ' ' (2.32) 

This definition is correct because each operator A 'e^ ( - tk + 1 ~ tk ^ HA , k — 0, 1, . . . , n, 
with to — and t n+ \ = [3, is of the trace class. The state (2.32) is obviously 
analytic in the domain 

< Re t 1 < Re t 2 < ... < Re t n < J3 

with boundary values at Re ti — 0, % — 0, ...,n, and it is uniformly bounded 
by the generalized Ginibre-Gruber inequality [23]: 

n n 

;=o i=o 

Here a^(B) is the C*-automorphism of C(TCa)- 

(4(B) = e- irAA Be tT " A , B e C{H K ) 
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and so the boundary value of the state p K satisfies the KMS conditions (see 
[2,32]). Moreover, the following formula is true: 



IV /ill/ 11/ 

1=0 Jn flA i=o 1=0 (2.33) 

where d[i A is the Gibbsian modification of the measure d/i A : 



rf/i A (0 A ) = Z^Ay'e-So^nM^d^A). (2.34) 

Following [7] we shall call this measure the Euclidean Gibbs Measure (EGM), 
which corresponds to our particular model (2.1)-(2.3) in this context. 

So, due to this construction the Theorem 2.1 for quantum Gibbs states 
(2.32)-(2.33) can be reformulated as follows: 

Theorem 2.2 For the system of quantum particles with Hamiltonian (2.1)- 
(2.3) there exists a sufficiently small mass m* such that for any < m < m* 
the weak limit of measures (2. 34) 

lim u A = a, 

exists and \x is a unique translation-invariant Gibbs measure on (£l* Z v, 
^g ■z v ) ^ n th e se t °f tempered Gibbs measures, for all temperatures T > 0. 

Remark 2.2 To prove the existence of the limit measure we are going to ap- 
ply a cluster expansions procedure to the perturbed measure /i A (2. 34), which 
is constructed by the Hamiltonian (2.1)- (2.3) with p.b.c.(2.4). But, in fact, 
one can prove the convergence (in the thermodynamic limit) of cluster expan- 
sions for some non-empty set of boundary configurations of so-called tempered 
type (see (3.34) of the next Section), and so to prove the existence of the limit 
measure in the class of measures supported on these tempered configurations 
(the set of tempered Gibbs measures) . 

To prove the uniqueness of the limit measure in this set we use one elegant 
criterium of uniqueness due to [6] (Theorem 3.1). 
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Remark 2.3 To prove the existence of the limit measure \x it is sufficient to 
prove the convergence of the averages (2.33) for local observables 

Ab{<I>b) = J[M<l>r i ) (2-35) 

j&B 

for any bounded subsets B C A and Tj G S*. 

Remark 2.4 Theorem 2.1 is equivalent to the Theorem 2.2 on the subalgebra 
constructed by (2. 29) -(2. 31). Then we use the fact that any operator A G 21 
may be approximated by a linear combination of operators of the form (2. 31) 

(see Theorem 2.2 of [32]). Hence the state p h (A) will also be approximated 
by (2.33). 

3 Cluster expansions and proof of main The- 
orems 

Remark 3.1 Before we construct the cluster expansions, which we need for 
the proof of Theorem 2.2 we recall the difference between the case s > 1 
in (1.3) and our case (1.5) with harmonic growth at infinity (see Section 1 
and discussion in [43,44])- Notice that for harmonic interaction (i.e. for 
polynomials of the second degree) the rescaling (2.6) implies, for s > 1, a 
week interaction between oscillators proportional to m( s-1 )/( s+1 ). But for the 
harmonic growth of the one-site potential and for harmonic interaction the 
rescaling (2.6) leads to smallness of anharmonic perturbation V A , see (2.10). 

Remark 3.2 Notice that the cluster expansions developed in [8,36,44] f or 
s > 1 are suited for small interaction between oscillators. In contrast to 
that the cluster expansions presented in the present paper are relevant for 
small one-site anharmonism and harmonic oscillators interaction. Therefore, 
it is suitable to arrange our reference measure in such a way it includes 
this interaction. Following the paper [2] it gives us the Gaussian measure 
representation for the state. 

The cluster expansions technique for the Gaussian integrals with non- 
diagonal covariance was developed in [29] (see also [30,49]). Then it was 
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applied in [20,22] for the proof of the Debye screening in dilute charged 
particle systems. This type of cluster expansions includes derivatives acting 
on the Gibbs factor exp[— V A ]. Since the coefficient b m is proportional to 
m 1 / 2 , and the coefficient 5 m is proportional to vrT 1 ! 2 (see (2.10)), we have no 
analyticity at m — 0. Therefore, the estimates are delicate in the small-mass 
domain. 

To prove the convergence of the measure dfi A in the thermodynamic limit 
it is sufficient to prove the existence of limit state for observables Ab (2.35). 

We would like to notice that similar to non-polynomial quantum field 
theory (see e.g. [30,49]) our cluster expansions correspond to a one-site per- 
turbation, and not to the interaction between sites. To start the construction 
of cluster expansions we define new configuration space. First we consider 
Sa for large (3 including f3 — oo. So for a fixed m we consider sufficiently 
low temperature such that critical behavior (if any) might be in the interval 
[0, p\. Then just for convenience we consider J3 G IN to define the partition 
of the interval [0, /3] into unit intervals. So, we define 



^ ■- 



-- {A f :=[f,f + l]|f = 0,l,...,/3-l, /3gw}, (3.1) 

Then we define 

t:=(j,r), je%», TeS $ , (3.2) 

and the box 

T:=A®S0CZ v ® M+, (3.3) 

where A and S& are defined in the previous section. With these notations we 
introduce 

^(M'):=G^,(r,r'). (3.4) 

For any j G 5Z U and Af G As we put 

A,,, :=(j,A f ), (3.5) 

and consider T as finite union of "rods" A^f, j G %", f = 0,l,...,/3— 1. 
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To construct the cluster expansions for the expectations of (2.35) we 
define (following [20]) a family of sets by the inductive procedure. Let for 
n — 1 

X 1 = Y 1 = (B,A B ), A B =UA„ (3-6) 

jeB 

where Aj coincides with Af such that Tj G Af . For n > 2 we put 

Y n =(j,A*),jeA\B, or A f e Ap \ A B , (3.7) 

and define 

X n := X n ^\jY n , X c n :=T\X n . (3.8) 

So, Y n for n > 2 are "rods" A^f, and X ra are their unions. 

Now define the sequence of new covariances (for < Sj < l,i = l,...,n): 



Gi (t, t'; (s) n ) = Gi (t, t')p(t - t'; (s) n ), (3.9) 



where 



n+l 



p(t,f;( S ) n ) = 5]ly m (t)ly T „(0 



m=l 



^ S iSi+1 • • • S m _! [lly ^ (t)lly m (t , ) + l Y Jt)l Yl (t')] . 

l<Km<n+l (3.10) 

Here ly(t) is indicator of the "rod" Y, and the "rod" Y n+ \ we identify with 
X c n . 

One can treat G^(t,t'\ (s) n ) as a kind of "interaction potential" in the 
space of "cluster" configurations A x A*. Then the interpolation parameters 
Si specify the intensity of this "interaction" between "rods" of Xi and of Xf. 
Putting s n = we get that there is no "interaction" between the "rods" of X n 
and "rods" of X^. Then let //jL s \ be the Gaussian measure with covariance 

<jq (£,£'; (s) n ). This is the kernel of a positive-definite operator in L 2 (T; lR d ) 
because it is a linear combination with positive coefficients of operators of the 
form 1_xGq1-x, where X is some of the unions of Y m , see [21], [30, Sec. 18.2] 

17 



and Section 4 for details. It is clear from (3.10) that A^>(i) n = A*o • F° r s « = 
one gets that the measure djj^., 3 ^ o becomes the product measure 



A,X n ,, A,X£ 



"/ i 0;(s) n _i,0 — "/ i 0;(s)„-i ^0^ 



(3.11) 



A,X n 



where the covariance of du r < ', ? is 



G£' A "(M'; (s)n-i,0) = (#(*,*'; (s)„_i,0) - l X c(t)Gft(*,OW*') 



(3.12) 



and the covariance of d/i .' Q n is 



T/j. + A 



G J ^(t,t') = t x ^t)G 1 (t,t')l x .(t'). 



(3.13) 



Notice, also, that (s" 



j i •--'"^-'■i 



ra-1 



n=\ 



and du n \ ? 



each step in the expansion can be obtained by app 
formula: 



A,X 



n=\ 



: d/i ' 1 , and that 



ying the Newton-Leibnitz 
(A B (-))^ = Z T X I dnl (<f>)A B (cj> B )e- V ^ + 



H /3,A 



Z T X \ dsi 



o d Sl Jn , A 



d&MMMe-*™ 



d4' Xl (<t>xMB(tB)e-^^?^ + Rt 



H 0,A 



where Z? '■= Z&(A) and 



Z T (X)= / d^(0)e 



T (rh\ e -n<t>x) 



H 0,A 



(3.14) 



(3.15) 



Using integration by parts formula to calculate the derivative w.r.t. param- 
eter Si (see, for example, [30, Sec. 9.1]) we obtain 



1* 



Rf = Z T l [ d Sl [ rf/i^ i (0)A^ B (0 B )e- f( ^ ) , (3.16) 

Jo Jh 0a 

where 

A, = X>i a) > ( 3 - 17 ) 



a=l 



A'f = ly f dT f dT'G£., k (T-T')4-p(U,T),(k,T'y,S 1 ) ,, S ' , , 

= : \L dt L MGlM ^ pltAsi) wS¥W) 

For definition of Ar^, (p,p' G IV) see (3.24). Then, taking into account that 



-v(4>t) _ P -v{<t>x 2 ) p -y(<t>x-) 



(3.19) 



with X 2 = X 1 UY 2 = FiUF 2 , and putting (3.17)-(3.19) in to (3.16) we repeat 
the procedure in each term of the sum over Y 2 to carry out the second step 
of expansion. After tit = \T\X\\ steps all "rods" in T will be exhausted 
and we get the identity: 



TIT 



(Ab(-)U = EE K$(Ab)FZ(Y), (3.20) 

ra=l YGT 
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where 



ICl(A B ) = !{ds) n - l f{ m (s) n - 2 )lT(Y, v ,A B ; (s)n-i), 

J (3.21) 



^,a (3.22) 

with 

n d 

A( V , Y) = ft A vm , A PtP , = J2 Ag" ( 3 - 23 ) 

/ = 2 Q = l 

where 77 := {(??(/)) OIL2} i s a sequence of couples with r)(l) < I, which we use 
to construct the trees estimates according the standard procedure, see e.g. 
[21], or [42] Ch.II.4. Here 

^W/'I/^ V'w-Wr (3 ' 24) 



ip " j p - 



/(??; (s)„- 2 ) = Y\ (Sr,(m)S v (m)+l ' "Sm-2), (3.25) 

2<m<n 

77(771) < m — 2, s = (si, . . . , s n _i), < Sj < 1. (3.26) 

The second sum in (3.20) is over all possible sets Y = {Y 2 , ..., Y n } in T. The 
product in (3.25) is equal to 1 if 77(772) = m — 1. Finally 

FZ(Y) = ^p-, Z T (X c n ) = Z T (T\X n ). (3.27) 

(For more details about construction of this type of expansions see, for ex- 
ample the survey [49].) 

To prove Theorem 2.2 we use the representation (3.20) and the following 
collection of lemmas. 
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Lemma 3.1 With the same assumptions as in Theorem 2.2 there exists a 
sufficiently small value of the mass m*, and constant C B which does not 
depend on T and on m, such that for all < m < m* 



J2 Kv(Ab) < C B e n -\ (3.28) 



where e = e(m) — ► as m — ► and 



lim K%(Ab)=K&{Ab). (3.29) 



Lemma 3.2 With the same assumptions as in Theorem 2.2 there exists a 
constant c such that 

F n T (F) < e^l*"" 1 ). (3.30) 

The constant c does not depend on A, ft, B, m, and one has the limit: 

\im v F?{Y) = Ft{Y). (3.31) 

We also need the following corollary of Lemmas 3.1 and 3.2. 

Corollary 3.1 Let B' C A and B fl B' = and let the sum over Y be 
restricted so that X n fl B ^ 0. Then 

oo 

J2 E \K Y (A B )\FT(Y)<Ci B e-^ dis % B ' B '\ (3.32) 

n=l Y,X n nB'^9 

where a is the constant of the one-site harmonic potential. 

Now, using the standard "doubling measure" trick by J. Ginibre [28] (or 
[30, Sec. 18.3]) we obtain from Corollary 3.1 the exponential clustering of 
correlations: 

Lemma 3.3 With the same assumptions as in Theorem 2.2 and for any 
measurable functions Ab 1 ,Ab 2 w.r.t. a -algebras E a Bi andT,p B2 respectively, 
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where B\, B2 C A and B\ D B2 — 0, there exists a sufficiently small value of 
m*, such that for < m < m* and any temperature one has the estimate 



(A Bl (-)A B2 (-)}^ - (A Bl (-))^(A B2 (-)}^ < 

< C 1 (^ Bl )C 2 (^ B2 )e-^* 5t(Bl ' B2) . (3.33) 

uniformly in A. 

The proofs of Lemmas 3.1-3.2 and Corollary 3.1 are presented in the next 
Section. They imply the proof of our Theorem 2.2 for the case of p.b.c. 

To prove the uniqueness of the limit measure in the Theorem 2.2 for 
non-periodic boundary conditions, we keep the DLR language (see [25]) and 
consider EGM, with some general boundary conditions £, from a class studied 
already in [6,7]. To this end we define the set of tempered configurations: 

fy : = U % = U U e %*" I E e ' M H6Hl"M < °° , 

p>o p>o I ie% v ) (3.34) 

Remark 3.3 For tempered configurations Q^ defined in [6] one has: Q f ~ C 

VL l ~, since we include also configurations that may have an exponential growth 
at infinity. 

First we define the harmonic measure with zero boundary conditions 
Hq(- I 0) in the same way as it was done in (2.23) for /!q (• | p). Then for any 
£ G VL P ~ with p < y/a, where a is the one-site harmonic constant from (2.8), 
we define the perturbed measure with non-zero boundary conditions by 



(3.35) 
where 
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<9A 

£(...):= £ (-)■ ( 3 - 36 ) 

(l,l') ieA,i'eA c ,\i-i'\=i 

We also need the following auxiliary measure /i A (- | y\) on Wa A , which 
depends on some fixed trajectories y\ G fi^ A : 

d/2 A (x A | y A ) := %A | y^e' ^ dT ^ x ^ y ^d^{x K | 0), 

(3.37) 

where (see (2.10)): 



V(x a (t) I |/ A (r)) =b m J2 [e-^'M+wM) 2 / 4 + e -M*»-y») 2 /4" _ 

j€A (3.38) 

It is clear that with the measure ji (• | ?/ A ) one can construct the same cluster 
expansion as in (3.23), but with additional dependence on configurations y\: 



nj- 

(*M-)WbA) = EE *?(•** I y^ F n(Y I Va), (3.39) 

n=l YcT 

In the next Section we prove the following lemma: 

Lemma 3.4 Let m < m* as in Lemmas 3.1-3.3. Then for the measure 
fx A (- | y A ), defined by (3. 37)- (3. 38), the estimates (3. 30)- (3. 33) are true with 
the same constants which do not depend on configurations y A . 

Proof of Theorem 2.2. 

Existence. First we prove the existence of measures with p.b.c. and zero 
boundary conditions As we mentioned in Remark 2.1 the covariance C A . fe (r— 
t') converges pointwise to the translation invariant function Go-\j~k\{T — t') 
and the corresponding measures du,Q and da^fv \ are defined on Ea. So, 

^)V 5 /(n — 1) P 
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using Lemmas 3.1 and 3.2 we can rewrite the expansion (3.20) in the ther- 
modynamic limit: 

(^(•))m = E E i4(Ab)f£(X), (3.40) 

with the same definitions of I n and A(rj,y) (see (3.22), (3.23)), but with G 

and ^fc)n-i instead of G o and A*o;(S,-r 

Hence, there exists probability measures //# indexed by 5 C -ZJ 1 ' such 

that 

(^b) = / A B ((f>)dp B ((j)). 



JH §,A 

All these measures are consistent and consequently by the Kolmogorov The- 
orem one can construct a unique limit measure p on the a-algebra S^ z „ . 

D 

Remark 3.4 In the case of general boundary conditions ^ E Ci p - with p < 

y/a, the factors K,^, Ffi, and the constants c, Cb in the Lemmas 3.1-3.2 will 
depend on £. This means that from convergence of the cluster expansions we 
can obtain only existence of the limit measure, but not uniqueness. 



Uniqueness. To prove the uniqueness we use the following 

Proposition 3.1 (see Theorem 3.1 of [6]) 

Suppose that for every l G 21> V and r G [0,(3], for every increasing sequence 
of boxes A n C Zi v (A n — ► Zj" '), such that l G A n for all n, and for any two 
configurations C,,i] G VL P ~, p < \fa, one has the convergence: 



(0/o( r o)) M A "(-IO - (tkiTo))^^) -*■ (3.41) 

as n — > oo. Then there exist a unique limit measure \x which does not depend 
on configurations £,77. 

This statement is a quantum version of some known previous results, 
namely, the monotonicity arguments [41] for the measures on the configura- 
tions in the Classical Statistical Mechanics and a priori estimates for EGM 

24 



due to [10,11]. These arguments are true for the wide class of anharmonic 
potentials with ferromagnetic type of interaction. We use this statement as 
follows. Consider 



(<^o( r o)),uAn(.|£) - (0i„( r o))^A„ ( . 



M 



1 d 

ds-^(<f)i ( r o)}^(-\v+sO> C = t-V- ( 3 - 42 ) 



After differentiation we obtain: 



I(<^o( t o))mM-|£) - (<^o( r o)) M M-|»?)l < 



j J^ rP 

SU P o E / dri ^«o ( r ' r o I ^ + <)0'(t), (3.43) 

.■e[o,i] ^ „„Jo 



<M'> 
where 



(3.44) 

It is easy to see that the convergence (3.41) follows from the 

Lemma 3.5 For m < m* there are constants Ci,i = 1,2 independent of 
t, T',A,r], such that 



Kfr(i)<Cie-^ l - l '\, (3.45) 

where 

K*(l) = K£,(t, t> I r,), and K^(2) = [ dr'K^r, r' \ r,). 

Jo (3.46) 
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Finally, applying (similar to [6]) the Cauchy-Schwarz inequality in (3.43) and 
using (3.45), we get (3.41). □ 

Proof of Theorem 2.1. 

As we mentioned in Section 2, our Theorem 2.2 is, in fact, a reformulation 
of Theorem 2.1 on the algebra which is constructed by (2.29)-(2.33). The 
uniqueness of the limit Gibbs state on the corresponding algebra follows from 
(3.43)-(3.46) by virtue of Proposition 3.1 for the (fi^ v , Es %„) on the set of 
tempered configurations (3.34). The convergence of cluster expansions gives 
also analyticity of the limit state in the circle \m — ra*/2| < m*/2 of the 
complex plane m. The end of the proof is based on the same arguments as 
in Theorem 4.1 of [2]. □ 

4 Convergence of the cluster expansions 

To prove the convergence of the cluster expansions (3.20) for A /* Z u we 
should first give demonstrations of Lemmas 3.1 and 3.2. 
Proof of Lemma 3.1. To simplify the proof we omit indexes a = 1, ..., d 
and we take into account that the number of derivatives for each order of n 
should increase as d n ~ x . To write down the general term for any fixed tree rj 
(see (3.23)) let us introduce (following [16,17]) the numbers d v (k), k — 1, ..., n 
which are the characteristics of tree 77, i.e. d v (k) is a number of lines, coming 
into the vertex k from vertices m (m > k). It is clear that n-y = 0^(1), 
n k — d v (k) + 1 for k = 2, ...,n, are exactly the number of lines coming 
from vertices k — 1, ...,n and these numbers coincides with the numbers of 
derivatives supported in Y±, Y 2 , ..., Y n . Hence, 

n n 

^d rj (k)=n-l and ^n fc = 2(n-l). (4.1) 

fc=i fc=i 

Let us rewrite the operator A (77, Y) in the following form 



n n k r n 

iT/Jvk) +{n k ) 



k=\ 1=1 jYk k=2 



k ! X 



n 



26 



5 nk 

(4.2) 



Here we put 



fc-i 



i]2 := 1, r}k := 1 + ^^(0,^), & > 3, 



i=2 



where 8 n {i),r){k) is the symbol of Kronecker. In the following we also put 

A n k) . 

t k ■- t k . 

The main technical estimate is the following 

Proposition 4.1 Let D 1lk be parts of operator A(r),Y) and factor e~ v ^ Xn ^ 
corresponding to the k-th vertex on the tree rj: 

^n k — ^'n k \P r) ^y*li •••i*d JI {k)) = 



^o lS(fc) ~ t k J^o \h - t ii ' • ' G o Kh - t d v (k)) x 

{Y k ) n k 

x -^ r -e- v ^\ (4.3) 

Then for b m < 1 the following estimate is true: 

D nk (t^t u ...,t Mk) )<^b m ^ 2 n k \^ 

x sup / ^G^(t^-4- ) )G- T (.-t 1 )---Q T (--^ (fe) )e- 5 ^) 2 / 4 , 

(*)d v (k)JY k (44) 



"v 



here supremum is taken over all terms in the product, where one can choose 
G*q (• — ti) equal either to G^{tk — U), or f Y dt'G^it' — ti). 

Proof. For simplicity we demonstrate this proposition for the case when in- 
stead of functional derivatives we have the ordinary derivatives, since exactly 
the same line of reasoning gives the proof for functional derivatives. 

Let X = —b m e~ 5mX / 2 denote the "scalar prototype" of our anharmonic 
potential in (4.3), supported in Y k . Then it easy to estimate by induction 
the value of n-th derivatives of X: 

\X^\<b m 2 n 5^ 2 V^.e- Smx2/4 . (4.5) 
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This inequality can be obtained by induction using the recursion relation: 

I n (x) = -<J ro x/ n _i(x) -{n- l)S m I n _ 2 (x), I n (x) = e 5 - x2 / 4 X^(x), 

(4.6) 

and the inequality: 

\ x \k e -5 m x>/4 < 2^5-^^., . 

Now using the Leibnitz formula for the product X ( - 1 'e x , where the first deriva- 
tive X^ 1 ' = b m 5 m xe~ SmX < 2 , and the estimate (4.5), again by induction one 
obtains: 

e x \ <2 n b m 5^ 2 n\e- 5 ™ x2/A . (4.7) 



dx 
We also used here the fact that e x < 1. 

Remark 4.1 The inequality (4-7) explains the provenance of the coefficient 
in (4-4)> since by the same way this coefficient appears when one passes to 
functional derivatives (as in (4-3)), instead of the ordinary ones. 

Remark 4.2 Notice that we also should consider the terms in which some 
derivatives act on the factor Ab- But, since As is supported only in Y\ = 
X\, then even in the case, when for a given tree r] one has all r)(k) = 1, 
for k — 2, ...,n, and (n — 1) derivatives act only on As, the other (n — 1) 
derivatives in A(r),Y) act on the factor exp[— V a ((f>x n )]- Hence, we obtain 
the powers of the small parameter b m , that we need. 

This finishes the proof of Proposition. □ 

As a result we can write for /J the following estimate: 

n-\ 

I^Y, v ,AB;(s) n ^)<C(A B )(8db m 5 m y^ 1 l[d 11 (ky.x (4.8) 

fc=i 

n „ „ n 

x JJ (sup / dt k Gl(t' k ,t k )) I d^U^) JJe-*-^) 2 /*, t k e Y k . 

The next step is to estimate the Gaussian integral in (4.8). 
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Proposition 4.2 Let 

» n 

In = / d^f s) l x (0xj II e-^^l\ t k eY k , n> 2, 

Jn 0,A fc =i (4.9) 

then there exists a constant c, which depends on parameters d, J, a such that 

AY i / iJa \d\x n \/4 
/„<e c ^ ( V ) . (4.10) 

Proof. Write the product in (4.9) as new a Gaussian integral w.r.t. the 
measure drf(z) = l/{2Tt) dn ' 2 exp[-z 2 /2], z G R dn : 



a, a 



Interchanging integrals one obtains 



H , , JM an (4-11) 



A,X„ 



f n — d'y(z)e~ 5m / 4 ^- k '= lZkZk ' Go ' "(**>**' ;( s )«-i>°) (4-12) 



where GV "is defined by (3.12). By construction (3.9)-(3.10) the function 
GV "(^fc, ife'! (s)fi-i) 0) is a convex combination of "diagonalized" covariances 
(see [21], or [30, Sec.18.2]): 



Gi> Xn (t k , t k r, ( s ) n _ 1; o) = J2 x i°o\ E A > = *' ( 4 - 13 ) 



where coefficients Aj are monomials of the form s^ • • • s& (1 — SjJ • • • (1 — sj ,' 
and 



G? = E^w04(«). |JZ fc (*)=X n for all i. 

k k (4.14) 
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Substituting (4.13) into (4.12), applying then the Holder inequality, and re- 

\,x 

gets: 



turning back to the measure d\i ' w n , defined by the covariance (4.14), one 



/„ < n ( / d ^ e 

v Jm dn 



S m /^l k , =1 Z k Z k ,G^ Xn (t k ,t k ,)^ X ' 



n ( / ^(wn e_w(,i)l f = : n^ ( 4 - is ) 

i ^ W /3,A ° fc=l ' i 

By virtue of (4.14) we have: 

*i=n(/ <« fc(i) (^) n ^ m ^ )2/4 )-n^(^«)- 

* K \a ° yezlw 7 V (4.16) 

By definition (2.23) one gets for A = Z k (i) and F(0 Zfc(i) ) = riy eZfcW e-" 5 ™^) 2 / 4 : 

J k (Z k ( i ))=C (Z k ({);J)- 1 [ faz^e-i 1 ™*"®-**-*' ]J e~^^l\ 

JC1 0A Y&Z k (i) (4.17) 

Here the measure du>z k (i) is defined by l.h.s. of (2.23), with J = in definition 
of operator B A (2.15) and with A k (i) := Z k {i) fl A, S k (i) := Z k (i) fl S&: 

te fc(0 (...) = / dx A / < ; , A (<Me-i /„^(«A WlWA (r)HilVB A0 _ 

^, A • / ^ |A| • / n i8 ,A 

We put also : 



f dt(u(t) - ujiQf := J] /" 
fcW {l,l')cA k (i) Js 



dr(ui(T) -uj v (t)) 2 . 
s k (i) (4.i 8 ) 
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The normalization factor Cq := Co(Zk(i); J) is the integral (4.17) calculated 
for 5 m = 0. By the Cauchy-Schwarz inequality we get from integral (4.17) 
that 

uv KK h ' Yez k (i) Jil H (4.19) 

The factors Co for J and 1J can be calculated explicitly using (2.24), (2.18). 
The integral in the product (4.19) can be estimated (for |y| = 1) by 



du Y e' 5mUj{tY)2/2 = C (Y,0) / d/iL(0y)e- 5 "^W 2(Ty) 



^/3 /3,A 



J=0 



5 m e ^ + l\-"/2 mns , lvl / 2JE \d/2 



Co(y,0)(l + -%=.C^i)- fl/ < e const„|/ 



2v/a e $V^-lJ V2 v /a + 5 m / 4 2 q) 

This finishes the proof of the Proposition. □ 

Now we return to the proof of Lemma 3.1. Since 



/ dtGl{t' - t) = f dtGl{t' - t) 
Jy ' Jt 



Y Jy Jt jeA J ° 



=:J2 D o(i-3)< ^ D 0t {j) =: C G , (4.21) 

then taking into account dependence of b m and 5 m on m (see (2.10)) and 
definition (3.21), we obtain by summing over trees r\ the estimate: 



\^{A B )\ < C\A B )(l6b^C G m d / 8 ) 

^ p=l - 70 



n— 1 
X 
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The last step of the proof is due to the Battle-Federbush inequality [16,17]: 

n = \i]\. 
Jo 

This yields (3.28) with 



n „! 

EII d » ! / W'Vfa; (s)„- 2 ) < 4", 

r? p=l ^° 



e(m) = 646v^C G m d/8 . 

Now taking into account (3.30), we obtain the value of the small-mass thresh- 
old m*: 

m+ = (64b^C G e c y 8/d . (4.22) 

□ 
Proof of Lemma 3.2. The proof is standard. One can use general theory of 
polymer-type expansions [42] (or, equivalently, the Kirkwood-Salsburg type 
technique [30, Sec. 18. 5]), see, also, [49]. □ 

Proof of Corollary 3.1. If we restrict the sum in the cluster expansion 
(3.20)-(3.27) in such a way that, for example, Y<i C B', then all terms in this 
expansion are proportional to Gfo(ti,t2), where t\ G B, t 2 G B', then (3.32) 
follows from the asymptotic behaviour of Gou_ji(t) in the variable \i — j\. 
Note, also, that in the case of p.b.c. (see e.g. [30]) 

G 0;\i-j\( T I P) < G 0;\i-j\( T I °)> 

and it has exponential decay of Go-\i-j\(r | 0) in \i — j\, which is uniform in 
r. D 

Proof of Lemma 3.4. To prove this lemma one has to check that Proposi- 
tions 4.1 and 4.2 are true for the auxiliary potential (3.38). It is easy to verify 
that (4.5)-(4.7) are true for X± = X±(x | y) = -b m e- &m{x±y ^/ A with S m /2 
instead of 5 m . So, the estimate (4.4) is the same but with \e~ 5m{ - x ^ +y( ^ k ^ / 8 + 
e -5 m (x(t k )-y(t k )) /8] i ns t ea d f e -Sm,<P(t k ) /4_ Then in the proof of Proposition 

4.2 we have to modify the integral (4.20) to 



d[i^(x Y )e-^ x ^)±yrM) 2 , (4.23) 



H 0,A 



J=0 
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which gives 



/ s , \-<V 2 ^^ 

C (y,0)fl + — G#(0)J e «(i+* m G^(o)/4) 



j=o 



Since the last exponent is less than one, the dependence on the configuration 
y disappears, and this finishes the proof. □ 

Proof of Lemma 3.5. We follow the line of reasoning in the proof of 
Theorem 3.2 [6]. Let d/#(0, <p' | 0) := d^{(j> | 0) ® df4((fJ | 0). Then we can 
rewrite (3.46) as 



(4.24) 
After change of variables: 

by definition of the measures it is clear that 

dp,£(<l>, 0' | 0) = dfio(x, y | 0). (4.25) 

Then, using definition (3.44) we get: 

Kfr(r, r'\r]) = Z(A | r])- 2 [ [ x^x^T^e-v^^ x 



Wa . JHa 



\9,A""V»,A (4.26) 

xei^')^ drvi{T)r}l ' {T) d~^{x,x' | 0). 



Define the measure 
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~^\dx K | y K ) = Z(y K )- l e-^^)d^(x K | 0). (4.27) 

Then 
K a (t,t' I 77) = Z(A | V y 2 [ drt(y A I 0)Z(y A )ei^%Io d ^>h'(r) x 



H /3,A 



x(x / (r)^(r , )) /i A ( .| s/A ). (4.28) 

Now the proof follows from Lemma 3.4, or more explicitly from (3.33) with 
Abx = %i and Ab 2 — xi>. □ 



5 External field analyticity of the Gibbs state 

In the previous sections we have proved uniqueness of the equilibrium state 
for our model (in the DLR-sense) in the set of tempered EGM, see (3.34). 
In this section we consider the model (2.1), (2.2) in the presence of external 
field h: 



H A (h) = H A + V K + ^h- qj) heR d . (5.1) 

jeA 

Here Hq and V A are defined by (2.1)-(2.3). After rescaling (2.6) one gets: 

H A (h) = H A + V A + J2h-Xj, h = m~ x l%. (5.2) 

jeA 

Then the following statement is true: 

Theorem 5.1 For the system of quantum particles with Hamiltonian (5.1) 
there is m* 7 such that for any < m < m*, T > and \h^\ < h (m*), a = 
1, ..., d, the limit 

exists. The limit p@' h (A) gives a unique state on the algebra 2lo, which can be 
continuously extended to the algebra 21. Moreover, the functions p^ ,h (A) , A e 
21 are analytic in the domain {h^ E C d \ \h^\ < ho(m*)}. 

34 



Proof. The proof is close to that of Theorem 2.1. The only remark one has to 
do is that in expansions (3.20)-(3.27) the factor exp[— V(<f>x n )] exp[/i f x <p(t)dt] 

appears instead exp[— V(<f)x n )} and the operator A(i],Y) in (3.22) acts on 
the product Ab{4>b) exp[— V(<f>x n )] exp[/i m f x <fi(t)dt]. When the derivatives 
from A (77, Y) act on the factor exp[h j x <f)(t)dt], the convergence is guaran- 
teed by the small factors h = m^^h from {h £ C d \ \h^\ < h^m^^a = 
1, ..., d}. In principle this implies a new value of the light-mass threshold m^: 



ml = mm{m*,m*(hC G e c+l y 4 }, (5.3) 

where m* is defined in (4.22), Cq is from the estimate of the product of 
covariances in A(r/, Y), e c is from the estimate of F^(Y), and additional e is 
from the inequality (cf. [21]): 



Y, [ (dsT-'fim (s)n-2) < e n , 



n = \r]\ 



But by virtue of (5.3) there is h (m*) such that m^ = m* for \h^\ < ho^m*), 
which finishes the proof. □ 

As a consequence of the above theorem we obtain the following result (cf. 

[43]): 

Theorem 5.2 Let A„ be the sequence A ± C A 2 C ... C A n ..., U^° =1 A„ = W , 
and let h = e\\h\\, where e is a unit vector in M. d . By 

a{0) := lim lim pf (-L £ Qj • e) (5.4) 

h— >0n— >oc \ A„ * — ' / 

we define the order parameter for the model (5.1). Then there exists mass 
m* > such that for any < m < m* and for all temperatures, including 
(3 = 00 ; the order parameter is trivial for any tempered boundary condition 
defined in (3.34), i- e - : 



a((3) = 0. (5.5) 
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Proof. The proof follows directly from the Theorem 5.1 and cluster expansion 
for the state 



|An| jeA„ 
Indeed, if derivatives from A(r], Y) act on the factor 



exp < / h ■ <p{t)dt 

then after thermodynamic limit the corresponding terms are proportional to 
h = m~ l l i h. If derivatives act on the factor <^exp[— V((f>x„)], then in the 
limit we obtain a Gaussian integration with zero mean and with odd numbers 
of fields. So, taking into account that (f' h is analytic in h (Theorem 5.1), in 
the limit h = we get zero value (5.5) for the order parameter. □ 

6 Convergence of cluster expansions and Gibbs 
state uniqueness for high temperatures 

The problem of the existence and uniqueness of the quantum state (EGM) at 
high temperatures for the model (2.1)-(2.3) can be solved in the framework 
of a general approach, developed e.g. in [13,14]. But there exist an indepen- 
dent question of construction of high temperature expansions for quantum 
state (2.5) and of the proof of its convergence in thermodynamic limit. This 
question is related to a specificity of the one-site potential (2.3) behaviour 
at infinity. For the quantum systems with unbounded spins the correspond- 
ing expansions were considered in [35,36,45,46,53]. All these results were 
obtained for the models with the one-site polynomial anharmonic potentials 
of the form (1.3). Then the convergence is insured by the small parameter 
p(s-i)/2s_ gj nce j n (2.1)-(2.3) we have s — 1, there is no such parameter 
in our case. We recall that 2s is degree of one-site anharmonic polynomial 
interaction in (1.3). 

In the present section we show that approach we developed in Sections 
3-5 (the choice of reference measure and the type of cluster expansion) can 
be adapted to obtain convergent cluster expansions in powers of f3 > 0, and 
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to prove analyticity of the state in domain \(3 — (3*/2\ < /3*/2 of the complex 
values of (3 around the point /3*/2 > 0. 

The construction of cluster expansions is absolutely the same as above, 
but instead of partition of the interval [0, /3] = Sp into unit intervals we 
consider it as one unit with variable r running in Sp. In another words, 
instead of "rods" A^f (see (3.5)) we consider for every site j G A we consider 
"rods" Aj t p := (j, Sp). Then the measure d[/,Q., s \ is constructed in the same 
way as in Section 3 with help of the covariance G^ k {r,T'\ (s) n ), see (3.12), 
where X n = and T = A . Then following the line of reasoning of Section 3 
we obtain: 



"A 



(^(•)>,a = EE ^(^b)F^(Y), (6.1) 

n=l YcA 

with the same definitions as in (3.24)-(3.30), and with only difference that 

, . rP rP x2 

A< a >=/ dr drG^y y/ (r,r')— p: p , (6.2) 

P ' P Jo Jo ' P ' V S<$(T)6<f)M(T') 

and Y p , Y p > are the sites of the lattice %> v fl A. The main result then is the 
following 

Theorem 6.1 For the quantum system with Hamiltonian (2.1)-(2.3) there 
exists a sufficiently small (independent of mass m) value of inverse temper- 
ature (3* such that for all < (3 < [3* the weak limit of the measures 

lim a = u, 

exists and [i is a unique translation-invariant Gibbs measure on (^lp^ v i ^p,z v ) 
in the set of tempered measures. Here /x A as in Theorem 2.2. 

Proof. It is very close to the proof of Theorem 2.2. The only difference is 
that in the estimate (4.20) one should take 

( 1 *m e^+l \-*/2 < / 20y/Z y/* < 2(3^ d/2 

\ 2^E e^-lJ -\2$y/E+sJ ~ 6 ( 63 ) 
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This gives for the small parameter of cluster expansion the value 

e = 64b^C G (3 d/2 . (6.4) 

Therefore, 

ft = (mV^C G e c )- 2/d , (6-5) 

□ 

7 Conclusions 

In the present paper we construct the Gibbs states for quantum crystal with 
nonpolynomial anharmonic potential (1.5) and with nearest neighbour har- 
monic interaction in the light-mass regime. These states are constructed on 
the algebra of quasi-local observables described in Section 2. Our analysis 
is based on a reformulation of the quantum Gibbs states problem into the 
study of Euclidean Gibbs Measures (EGM) [2,7,31,33]. 

The main technical tool we use to get the main result is the cluster ex- 
pansions for EGM. They allow us (Theorem 2.1) to prove the existence and 
uniqueness of translation-invariant EGM in the light-mass domain, for all 
temperatures T, including T — {j3 — oo). 

To prove the existence of a unique translation-invariant measure we start 
with periodic boundary conditions (Section 2). Then we prove convergence 
of the corresponding cluster expansions and the cluster property of the EGM 
in the light-mass domain for all temperatures T. This gives the uniqueness 
in the class of translation-invariant measures for small mass m < m* and all 
T > 0. Next we establish that convergence of the EGM cluster expansions 
and their cluster properties hold also for zero boundary conditions. This 
again implies uniqueness in this case. Finally, following the arguments of 
[6,11] we prove the uniqueness of EGM on the set of tempered configurations 
(Section 3). 
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